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ABSTRACT 


Overbars: 


A rapid technique for solving.the blade-to-blade 
turbomachinery flow problem has been developed. Ap- 
proximate governing flow equations, which include the 
effects of compressibility, radius change, rotation, 
and variable stream sheet thickness are solved using 
an Inproved panel method. The development and solu- 
tion of these equations are described. S^le calcu- 
lations are presented to illustrate the method's ca- 
pabilities and accuracy. 
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Integral equation solution singularity 
Influence coefficients 
stream sheet thickness 
meridional blade-to-blade surface 
coordinate 

blade shape surface normal 

radius from turbomachine centerline to 

blade-to-blade surface 

blade row pitch or blade spacing in 

radians 

blade surface distance 

blade-to-blade and midchannel surfaces 

absolute velocity of flow 

relative velocity of flow 

disturbance relative velocity 

transformed meridional coordinate 

transformed tangential coorddinate 

flow angle on the blade-to-blade surface 

blade row circulation 

tangential blade-to-blade surface 

coordinate 

fluid density 

singularity strengths 

orientation of blade surface elements 

used in the Integral equation solution 

angular velocity 


Subscripts: 


c known or constant value 

est estimated value 

i,j,k summation indices for singularity strengths 

u,i,o,d upstream, inlet, outlet, and downstream 
station of the blade row 

n,t surface normal and tangential components 

onset undisturbed flow 

rb integral terms contribution to solution 

velocity 

p arbitrary point on blade-to-blade surface 

m,f on blade-to-blade surface of revolution 

x,y on transformed blade-to-blade surface 

* transformed quantity 


( ) vector quantity 


( ) mean quantity 

INTRODUCTION 

The design and study of turbomachinery rely 
heavily on the use inviscid Irrotatlonal compress- 
ible blade-to-blade flow solutions. As a result many 
methods and computer codes have been developed to 
solve this problem. The intent In presenting yet 
another of these blade-to-blade solutions is to pro- 
vide the turbomachinery designer with a tool which 
has a better conR)inat1on of solution speed, robust- 
ness, and versatility than previous methods. 

The procedure solves the inviscid irrotatlonal 
compressible blade-to-blade flow equations on a sur- 
face of revolution. One three dimensional effect is 
modeled by including a variable stream sheet thick- 
ness. The governing equations are linearized by ap- 
proximating the compressibility effects. These ap- 
proximate equations are solved using an Integral 
equation technique (panel method). 

The method does have limitations and is intended 
primarily for use in preliminary design studies where 
a variety of blade shapes and flow conditions are to 
be analyzed. By using this technique the designer 
can produce blade shapes for analysis by more accu- 
rate, but presumably longer running, flow codes. 

Specific features of the method are as follows. 
(1) Computational time for a typical case run on the 
IBM 3C33 computer at NASA Lewis Research Center is 
less than 4 CPU seconds. After an Initial solution, 
additional solutions for variations in the Inlet flow 
velocity, inlet flow angle, outlet flow angle, rota- 
tional speed, or a combination of thes'* effects can 
be made at an additional expense of less than 0.5 CPU 
seconds per case. (2) The method can be applied to 
all types of turbomachinery and flow conditions. 
However, the quality of the solution does vary with 
the particular problem being solved. The method 
loses accuracy e; the flow becomes transonic, will 
not predict shocked flows, and does not accurately 
determine when a blade row is choked. In addition 
the method also gives poorer solutions in the uncov- 
ered portions of the blade passage for lower solidity 
blade rows, (3) Since a direct linear equation sol- 
ver is used, a solution is always produced even if 
the blade shape is poorly described or extreme flow 
conditions are given, (4) The method uses an inte- 
gral equation solution, so that the flow Is only cal- 
culated on the blade surface. The flow can be cal- 
culated at points In the blade passage, but this is 
inefficient, (5) The solution does not require the 
use of a computational mesh or grid, so that almost 
any blade shape or body can be analyzed. 
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ANALYSIS 

The general blade-to-blade surface geosietry Is 
shown In figure and an r-z cut through the blade- 
to-blade surface Is shown In figure 2, The surface 
In figure 1 corresponds to the Sx surface described 
by Wu In reference [1]. The governing equations de- 
veloped by Wu for Inviscid, Irrotatlonal, compress- 
ible flow In the $x surface can be written a$» 


Continuity (pbV^) ♦ -|jj (rpbVj - 0 (1) 


Irrotationality -L- (rV^) - 0 (2) 


W W 

tan (a) m m 

m X 

Using this transformation the governing equations 
become 


Continuity ^ (W^) + (W^^) - 0 (8) 

Irrotationality a I ) a dr 

17 \7b/ - alVit/ “2flr ^ 

(9) 


where V. and Vjj, are absolute velocity compo- 
nents, p Is the fluid density, b Is the stream 
sheet thickness, and r Is the local radius of the 
axl symmetric blade-to-blade surface* 

In terms of flow relative to the blades the 
equations become 


Continuity (pbW^) + = 0 (3) 

Irrotationality -|^ (Hj - ^ (rW^) - 2ar ^ (4) 

where the relative components, and Wjj, are given by 
'^m == 

(5) 

V, = ♦ ar 


The radius, r, can be transformed out of the 
governing equation by using the same transformation 
as Kurmanov, Podvidz, and Stepanov [2], 


Expanding the Irrotationality equation (9) gives 

( 10 ) 

The effect of the transformation on the governing 
equations has been to move the non-linear terms from 
the continuity equation (8) to the irrotationality 
equation (10) and to have the radius expllcity appear 
In a single term. 

Equation (10) is linearized by replacing (pb) 
and Wy in the last three terms by estimated val- 
ues, which are taken to be functions of x onl;,. 

With these approximations the governing equations 
become 

Continuity (W^) |jj- (W^) = 0 (11) 

Irrotationality (W^^) - |^ (W^) = 

■ ^ * 2o(pb)g5^r ^ 


y • 9 i ^ 
or 



Equation solution 

Since equations (11) and (12) are linear, 
superposition of solutions may be used 


^ j a_ (6) 

ay “ ’ r ax * am 


and 


where r is a function m only. In the present 
work, the velocity part of this transformation is 
modiHed to the following 


Wx ® *^P^m * ^y * *"phW^ 


(7) 



(13) 



where W and W can be considered as components 
^c ^c 


of a flow which is indepe.ident of y , and w^ 
and Wy are the disturbance to this known flow due 
to the^blade presence. Substituting equations (13) 
into the governing equations gives 


Note that the transformation preserves the relative W * constant (14a) 

flow angle from the m-f surface to the x-y surface. ^c 
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W • constant + 


i*(p^)est 

J est 


/o‘>>est 

“rb “I “> est “<•'»> 


est 


- 2oJl 




The above equations, (15a) and (15b), correspond to 


- 2ij (ob)gj^ r dr 

(14b) 

tion technique of reference [3] may now be applied 
directly using slightly modified boundary conditions. 

h ° 

(14c) 

Boundary conditions 

The boundary conditions for the relative 
transformed flow problem are 

h ° 

(14d) 

1. On the body surface. 


In this study a solution to these equations was 
developed by using the integral equation approach of 
reference [3]. In the notation of this paper, the 
onset velocity reference [3] becomes. 


^onset 


■^onset 


“onset 


ref. 3 


ref. 3 


onset 


s1n(a) 


.(pb)est 

"> ,5, «(•") 


est 


- r dr 


and the distributed singularity expressions of refer- 
ence [3] are used to satisfy equations (14c) and 
(14d), The resulting equations for the velocity com- 
ponents tangential and normal to the blade surface can 
be written as. 

n 

“t = “onset «s(e -a) + sin(e ) ♦ 

P P i-1 


n-1 


m 


j«l k«l 


“n ■ -“onset V“^ ‘ “rb * E 


i-1 


®*p1 “1 


n-1 

>1 


^''pj 


m 


(15b) 


where 


¥.n ■ W^(s) « 0 


2. At the upstream boundary, 

“x|^ • “ “onset<=°®(“> * 


^ - W^s1n(.J - “onset®’''^*^ * \ 


3. At the downstream boundary, 

“x|j = “d"°s(«^) * “onsetPP*<“) * ^ 

“y|, ‘ V’"<“d> ■ “onset^’"<“^ * \ 

4. Circulation condition, 

r - S[(rV)y - (rV)^,] - ^rV.dl 
or in transformed variables 


■&I u -(^)d * “ - 


'•d') 


and 


r .(£ 1 (U + pb or^) * (dx^+ dy^) 

Ah 


7x1/2 


'7H 

5. Far stream limits on w, 
’‘u**d 


*x ■ “x ® ** * " 

u d 


-W as Xy.x,, -> * - 


(16a) 


(16b) 


(16c) 


(16d) 


(16e) 


(16f) 
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The above conditions* along with the conservation of 
angular momentuir (rV^ - constant) in areas 
outside the blade row and equations (1S)« used to 
develop a system of equations similar to equations 
(17)-(19) In reference [3]. This system of equations 
Is solved by upper and lower triangular matrix 
decomposition. 

Kean flow and stream sheet thickness estimation 
solution of the governing equations re- 
quires that the integral terms In e tatlon (14b) be 
known a priori. These terms can be approximated by 
using numerical Integration and a one dimensional 
estimate for the flow through the blade passage. 

Three techniques have been used to estimate the 
one dimensional flow in this study. The first is to 
use the blade shape and passage geometry to calculate 
approximate flow turning angles, flow areas* and 
stream sheet thicknesses. Thesf quantities are then 
uied with a one dimensional continuity equation and 
tfie isentropic relations to calculate the mean flow 
velocity and density. The second technique is to 
specify or assume a mean flow velocity distribution 
through the blade row and to estimate the stream 
sheet thickness from the passage geometry (see 
figure 2). Isentropic relations are used to 
calculate the density. The simplest mean flow 
estimate is a linear distribution from blade inlet to 
exit based on the design velocity diagram. More 
complex distributions can be made depending on the 
type of turbomachinery being analyzed and one's 
knowledge of the flow. The third method is to use a 
mid-channel stream surface flow solution, which 
corresponds to an $2 surface in Wu's notation [1], 
Such a solution can provide both the stream sheet 
thickness and mean velocity distributions. 

All three techniques were used successfully in 
this study. The quality of the overall solution pro- 
duced by using any one of the mean flow estimates was 
dependent on the particular blade geometry being ana- 
lyzed. For highly guided passages like those found 
in centrifugal machinery, the first technique pro- 
duced very good results. In lower solidity blading 
where the flow angle can deviate significantly from 
tht mean blade angle, the second method proved to 
give the best results. The use of a mid-channel so- 
lution also provided good results for a number of 
cases, but its success depends on how well the mid- 
channel solution models the flow turning through the 
blade row. In this work the MERIDL computer code [4] 
was used to calculate the mid-channel solutions, rnd 
did provide good solutions for most cases. For flows 
where the mean flow is supercritical and the flow 
turning is controlled by the congress i b1 lity of the 
flow, only the second and third techniques gave rea- 
sonable solutions. 

Method accuracy 

In considering the overall accuracy of this 
solution, two sources of error are discussed. First 
is the approximations made in the governing equa- 
tions. Second is the solution technique used in 
solving the equations. 

By looking at the terms in equation (10) which 
are neglected or approximated in the final solution, 
the applicability and accuracy of the solution can be 
assessed. The only flow for which the solution is 
exact, in the sense that no approximation has been 
made to the governing equations, is incompressible 
flow with constant ctream sheet thickness. With the 
introduction of a var. able stream sheet thickness, an 
approximation of the Wy term is required. This 
ter.i makes significant contributions to the solution 


when problems with large or rapid changes in stream 
sheet thickness are analyzed. The accuracy of the 
term Is controlled by the estimate of Wy. At far 
stream stations the flow is modeled as being uniform, 
and the estimate of Wy as being independent of y is 
consistent. However, near or inside the blade row 
where blade surface curvatures strongly influence the 
flow, errors in the solution will occur. Compress- 
ible flow calculations require the inclusion of the 
change in densi y terms in equation (10). The more 
compressible the flow the more significant these 
terms become. In the solution equations used in this 
study it is assumed that the density only varies in 
the meridional direction. Wood In reference [5] also 
makes the Si.me assumption, and it is currently used 
in the TSONIC code [6] to calculate transonic flows. 
An estimate of tangential variations In density could 
have been Included In the analysis, but this would 
require an approximate solution for the blade 
loading. Calculations indicate that the term Is not 
significant for subcritical flows. However, 
inaccuracy in the solution will appear when high 
speed flows with large blade loadings are being 
analyzed. Such flow occur in supercritical blade 
designs. The remaining term for change in density in 
the iroridiondl direction is approximated by using a 
mean flow estimate to arrive at values of the 
tangential velocity con^onent and tne density. This 
term tends to dominate the solution as the compress- 
ibility effects increase. If the mean flow estimate 
does not approximate the local flow well, errors In 
the solution will occur. 

Major sources of error in the solution technique 
are the integral equation solution and the numerical 
approximation of the integral terms in equations 
(15a) and (Ibb). The integral equation solution is 
second order accurate in terms of the elements used 
to approximate the blade shape. Its accuracy was 
discussed in reference [3], Numerical integration of 
the terms in equations (ISa) and (15b) was accoirv- 
plished using a simple trapezoidal method. Twenty- 
two integration points were used between the upstream 
and downstream boundaries of the flow field with 
twenty of the points placed between the blade row 
leading and trailing edges. The accuracy of this 
integration is consistent with the accuracy of the 
integral equation solution for a typical design cal- 
culation. 

Solution modific at ions 

The solution can be modified to produce better 
results for certain problems. One possible modifica- 
tion is to specify a fictitious mean flow distribu- 
tion or stream sheet thickness in order to account 
for viscous flow blockage effects. Another solution 
modification is to alter the transformation back to 
physical variables. The velocity transformation can 
be written as. 




where W is the physical velocity magnitude and 
W* is the transformed velocity magnitude. This 
can be modified to a form similar to the 
Stockman-Lieblein compresslbil Ity correction [7], 



( 17 ) 
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where Is the transformed mean velocity magnitude at 
the s«ne merdlonal location as W, This altered 
transformation models density variations In the 
tangential direction which were suppressed by the 
approximations to the governing equations. 

Additional solution modifications are possible. 
The problem with usln^ solution modifications is that 
they usually produce better solutions only for a par-> 
ticular type of turbomachinery flow. 

SAMPLE CALCULATIONS AND DISCUSSION 

The sample calculation problems were chosen to 
demonstrate the capabilities and accuracy of the 
method. The calculated flows range from Incompress- 
iDle to '•upercritlcal. Both centrifugal and axial 
turbomachi nery blade rows are analyzed. All flows 
are calculated without taking Into account viscous 
effects. 

Incompressible flow 

Two inconH)ressible flow cases are presented. In 
both cases a centrifugal compressor blade-to-blade 
geometry is considered. The transformed surface 1' 
shown in figure 3. The blade row is rotating at 
68,000 RPM. In the first case a constant stream 
sheet thickness is assumed. Therefore no 
approximations are present in the governing 
equations. The calculated results are shown in 
figure 4, along with results from the finite 
difference solution of Katsanis [6]. The calculated 
surface velocities from the two solutions are in good 
agreement, including the region of reverse relative 
flow on the blade pressure surface. The second case 
differs slightly Irom the first. A varible stream 
sheet thickness is assurr^d, so that the method is 
solving approximate governing equations. The mean 
flow is estimated using the blade passage geometry. 

The assumed stream sheet variation has a ratio of 
4.75 to 1 from the blade leading to trailing edge. 

The comparison of surface velocities with Katsanis 
[6] is shown in figure 5. Agreement between the two 
solutions is very good. 


Compressible flow, subcritical 

Three subcritical flow cases are presented. The 
density approximation will influence the solutions 
for these flows. The mean flows used in the 
solutions were determined from the blade passage 
geometry in the first and second cases, and was 
specified in the thiro case. 

The first calculation is for a radial diffuser. 
The critical velocity ratio at the blade inlet is 
0.8/ and at the outlet is 0.41, As shown in figure 
6, the comparison with Katsanis is excellent except 
at the trailing edge. The trailing edge velocity 
spikes calculated by the Katsanis solution are due to 
the rapid turning of the inviscid flow around the 
rounded trailing edge of the blade. The real flow 
would separate forming a viscous wake before such 
large velocity gradients could appear. In the pres- 
ent solution an extrapolation technique similar to 
Gostelow's [8] is used to eliminate these spikes. 

The second case is a radial inflow turbine de- 
sign by Klassen using the ^ERIDL and TSONIC computer 
codes [4], [6], The critical velocity ratios across 
the blade row are G.48 at the inlet and 0.88 at the 
exit. The flow incidence with the blade is 60 de- 
grees, Total turning of the flow is 100 degrees. 

The rotational speed of the rotor is 1(X),000 RPM. A 
comparison of the present method with the TSONIC code 


solution 1$ shown In floure 7. Some d1fferem:es In 
the solutions do occur In the uncovered portions of 
the blade row at the Inlet and outlet areas, but 
overall agreement Is good. 

The third subcritical case Is a high turnina low 
solidity turbine vane. The vane has an axial solid- 
ity of 0.63, and a design turning of 75 degrees. The 
Ideal exit critical velocity 1$ 0.896. Since the 
flow turning deviates significantly from the mean 
blade angle for this problem, It was necessary to 
specify the mean flow rather than calculating It 
using the blade passage geometry. A linear mean flow 
distribution was assumed based on the vane velocity 
diagrams. The etperimental data of Schwab [9] Is 
used for comparison in figure 8. The agreement Is 
good. 

Compressible flow, supercritical 

Two supercritical cases are presented. The 
density approximation will have more of an effect In 
these cases than It had on the previous solutions. 
Solution modification techniques were used In 
calculating these flows. 

The first case is a turbine vane which was de- 
signed by Detroit Diesel Allison and tested there by 
Huffman et al. [10]. The mean flow for this calcu- 
lation was determined by using the blade geometry 
from the inlet to the throat and a linear distribu- 
tion from the throat to the blade row outlet flow 
conditions. The comparison with experimental data Is 
shown in figure 9 for both a supercritical and sub- 
critical exit flow. The agreement for both cases is 
good. The calculated velocity is a little low on the 
suction surface near the throat and blade trailing 
edge. 

The final case is a supercritical compressor 
blade designed by Sanz using the hodograph method of 
Bauer, Garabedian, and Korn [11]. The critical ve- 
locity ratio is 0.74 at the blade inlet and 0.58 at 
the outlet. The flow is turned through 29 degrees by 
tiie cascade. The mean flow for the solution was cal- 
culated using the blade passage geometry. This case 
is a severe test of density approximation, since 
large blade loadings and high speed flows are pres- 
ent. Solution comparisons are shown in figure 10. 

The circle symbols show the calculated results when 
the basic method is applied. The calculation does 
poorly in the region where the flow is transonic. 

The calculation in this region can be improved by 
using equation (17) to modify the solution as shown 
by the square symbols in the figure. However, the 
subsonic portion of the solution becomes poorer as a 
result. 

As this case demonstrates, solution modification 
techniques can be used to match almost any known 
blade row flow. The difficulty with using solution 
modifications lies in their general application to 
the prediction of flows. Hopefully, speific modifi- 
cation techniques could ^e developed and applied to 
the prediction of flows in particular classes or 
types of turbomachinery. 


CONCLUSION 

The present method rapidly calculates surface 
velocities for a wide range of turbomachinery blade- 
to-blade flows. The method produces good results for 
subcritical flows. Supercritical flows can be calcu- 
lated, but the solution accuracy is poor. Solution 
modifications can be made to improve the supercriti- 
cal flow calculations, however, application of these 
modifications requires insight into the flow being 
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analyzed and the modifications can not be applied 

universally to all flows. 
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Figure L - Blade-to-blade surface of revolution. 



CENTER LINE AXIS-Z 

Figure?. - R-Z cut through blade-to-Made surface. 
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Figure 4 - incompressible flow - constant stream ihi 
thickness. 
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MERIDIONAL COORDINATE 

Figure 5 - IncompressIMe flow-variable stream sheet thickness. 
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Figure 9. 




